Quasi-particle characteristics of the one-dimensional polaron 



G. Ganbold * 

Bogoliuhov Laboratory of Theoretical Physics 

Joint Institute for Nuclear Research, Duhna, 14-1980, Russia 

(September 15, 1999) 

0\ ; Abstract 

a^ 

CJ ' The main quasi-particle characteristics of the one-dimensional polaron are 

^^ ' estimated within and beyond the most general Gaussian approximation at ar- 

bitrary electron-phonon coupling. We have derived explicitly the ground-state 
^^ \ energy and the effective mass in the weak- and strong-coupling regimes. For 

arbitrary coupling, the Gaussian leading-order term of the polaron self energy 
^ ' improves the corresponding Feynman estimate and represents the lowest up- 

^vq ' per bound available. We have calculated the next non-Gaussian corrections. 

cn . Taking into account systematically higher-order corrections does not perturb 

considerably the obtained results. 
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I. INTRODUCTION 

A conducting electron moving slowly in a polar crystal interacts with lattice vibrations 
and forms a quasi-particle, which is commonly called a polaron. Traditionally, the polaron 
problem has been considered in three dimensions Hi]]- 0- In recent years, polaronic effects 
have been observed in low- dimensional systems [0-0 that has stimulated many theoretical 



enhancements of the conventional theory by lowering the spatial dimension |10|- |12 



The polaron confined in one dimension (d = 1) is usually modeled by the following 



Hamiltonian [ p^ , p^ 



^ = y + E 4a. + /f E (4e-^''^ - a,e^'^) , (1.1) 

k k 

where p and r denote the momentum and position of the electron; /c, a^ and a\ are the wave- 
vector, annihilation and creation operators of a phonon; L is the length of lattice crystal 
and the interaction is characterized by the dimensionless coupling constant a > 0. Here we 
are using appropriate units, such that, m?, = ujlq = h = c = 1, where rrih is the electron 
bare mass, ulo denotes the constant frequency of the longitudinal optical phonons. 

Various methods have been devoted to investigate the properties of the system described 
by ( |1 ■ 1| ) , in particular, the ground-state energy (GSE) and the effective mass (EM). However, 
strict results for these quantities are known only in the limiting cases of the weak- {a -^ 0) 
and the strong-coupling (a — > cx)) regime [p!0|, p!T| , p!^ . 

Among the variety of approaches providing reasonable results for finite a, the Feynman 
path-integral approach |^ stands out for its flexibility and all-coupling nature. Following 
this formalism one can integrate out the phonon variables in ( |1 . 1|) and define the free energy 
of the polaron J?-/3(a) as follows: 







where f5 is the inverse temperature, So[r] represents the free part of the polaron action while 
the interaction is described by 



c^W[r] '^^Jjdtds- :^^^ S{\r{t) - r{s)\) . (1.3) 







Functional integration symbol §6r in (|1.2| ) indicates path integration performed over all 
closed paths r(0) = r(/3) = x followed by ordinary integration over x. Hereby, the free 
energy of the non-interacting phonons is already subtracted, so the standard normalization 
condition J-'/3{0) = is obeyed. 

Thus, the original many-body model ( |1 . 1| ) has been transformed into an effective one- 
particle problem. However, the nonlocality and ^-singularity arising in W[r] prevent any 
further analytic treatment for finite a. For higher dimensions {d > 1), another Coulomb- 
like singularity arises in ( p..3|) . Among the approximate methods, Feynman's all-coupling 
variational estimate M has been particularly successful (for d = 3) to interpolate smoothly 



between the weak- and the strong-couphng regime of the ground state {(3 — ^ oo) by intro- 
ducing the following retarded oscillator 



/3 

aWplr] = -Cjjdtdse"''^'-'^ (r(t) - r{s)f (1.4) 



with adjustable positive parameters C and uj. A combination of the Jensen-Peierls inequality 
with the variational principle leads to the Feynman upper bound to the exact GSE. A similar 
scheme has also been suggested by Feynman to estimate the EM. Further improvement of 
Feynman's method has been made (for d = 3) hj increasing the number of trial oscillators 



An essential generalization of the Feynman method has been proposed in |T^ and |T^ 



independently. In |jT^ the Feynman exponential C exp{—uj \t — s\) has been changed to 
an isotropic trial function f{t — s) > obeying the periodical condition f{t — [3) = f{t). 
Appropriate variational optimization has lead to a coupled integral equation for function 



f{t). By solving it numerically one obtains the best upper bound available for rf = 3 [17 
According to |T3] this result can be re-scaled into another spatial dimension. 



An improvement over the Feynman estimate in a nonvariational way has been made by 
applying the Gaussian-equivalent representation method |T^J2^. According to this method, 
the polaron action in arbitrary dimensions d > 1 can be identically transformed into an 
equivalent form 

SGER[r] = -(r, D-h)+ : ^sM : , (r, D-'r) = jj dtds r{t)D-\t, s)r(s) , (1.5) 



where the Green function of the differential operator D^^{t,s) should be derived from a 
constraint equation which ensures the normal-ordered form : 5*2 [r] : of the new interaction 
functional and the absence of any quadratic path configurations in it. Within this method 
the entire Gaussian contribution and the next non-Gaussian correction to the GSE have 



been calculated for arbitrary coupling in two and three dimensions ||21| , p2| . The Gaussian 
term for the GSE represents an upper bound and in particular, for d = 3 we observe exact 
coincidence with the AGL-S result. Appropriate estimate for the EM within and beyond 
the GER method has been performed in [^ . 

In the present paper we develop the most general Gaussian approximation for a wide class 
of path integrals and apply it to the specific case of one-dimensional polaron by estimating 
its main ground-state characteristics and improving the known results. 

II. POLARON GROUND-STATE PROPERTIES 

The main quantities of interest, characterizing the quasi-particle properties of the polaron 
are the GSE and EM considered at the zero-temperature limit f3 -^ oo. The GSE of the 
system is 

E{a) = lim J^«(a) . (2.1) 

/3— >oo 



The EM of the polaron may be defined in different ways. To define the GSE and EM si- 
multaneously we consider the polaron partition function projected at small fixed momentum 
p as follows: 



-dJ'pip^a). 



V^ 



r(/3)=x 



-f3Ep{a) 



(3^00 



(2.2) 



r(0)=0 



Normalization in (|2.2|) is chosen so that for large (3 we obtain Tp{p'^, 0) = /3p^/2. By changing 
the integration variable r(t) — >■ q{t) + r{t) with q{t) = xt/ (3 one goes to conventional closed 
paths starting and ending at zero and rewrites ( p.2|) as follows: 



-I3rp(p^,a). 



CXD 

lirp J 



V^ 



(2.3) 



r(/3)=0 



{a,x)= y'5re-^«M+°^['-+'?l] = y't/ao[r]e'^^["+'?l, Jdao[r]-l = l. (2.4) 



Zpiia^x 



^/3l«, 



r(0)=0 



Note, going to the Fourier transform for (5-function we can rewrite the polaron self- 
interaction in (11. 31) as follows: 



/3 



W\r] 



2^2 



7r 



dtdse-^'-'^ /^^e^fc(Kt)-Ks))<^J:' dn,{t,s,k) e^'^'^^'^K 



(2.5) 



Since the polaron action is translationally invariant, Ep{a) is a continuous function of p 
and one can expand it around small momentum as follows: Ep{a) = E{a) + p'^/2m*{a) + 
O(p^). Therefore, the GSE and EM of the polaron are 



E{a)=— lim JFg(0, a) , 

/3— »oo 

m(a)=\ lim — -J^^(p^a) 



(2.6) 



p=0/ 



III. GENERALIZED GAUSSIAN APPROXIMATION 



As is known, no exact evaluation of (|2.4|) is available yet. Our strategy is to extract 
exactly the most general Gaussian contribution out of Zp{a,x). Then, the remaining non- 
Gaussian correction may be systematically estimated to improve the main approximation 
and to control the accuracy. 

i. First, we demonstrate the basis of our method by evaluating Z^(a,0). This simplifi- 
cation does not harm the GSE and the necessary extension to x 7^ required for the EM 
will be given later. In principle, the interaction functional W[r] may be more general than 
that given in (|1.3|) and (pT5|). So, the aim is to find an optimal representation for 



Z^(a) = (e"^W)o, (•)o = /rfao[r]., (l)o = l. (3.1) 



For a -C 1, the Gaussian contribution in ( |3.1|) is mainly represented by the functional mea- 
sure (io"o[r] and the influence of the quadratic part entering W[r] is insignificant. However, as 
a increases, the nontrivial Gaussian part of W[r] plays considerable role and this drastically 
diminishes the efficiency of the original representation ( ^.1[ ). Besides, the first cumulant 
(iy[r])o is nonzero. 

To describe the system more efficiently at finite and large a we go to another functional 
representation based on the most general Gaussian measure da[r] as follows 



(1) = 1, (3.2) 



A specific restriction imposed on the Green function D{t,s) of the operator D ^(t,s) will 
be given later. Obviously, for /5 — > oo we have 

(r(t) r{s)) = D{t - s) , (e^'^^^*'^)) = e-^'^^*"^!) , F{t) = D{0) - D{t) . (3.3) 

Functional averaging schemes (|3.1| ) and ( p.2|) are related to each other as follows: 



r(/3)=0 


= (•), 


r(0)=0 





(e^H.) = (e^H)(.)^ = e-^o(.)o, (3.4) 

where 



T[r] '=' \ (r, (D- - D,^)r) , C, '^ - In ^^^ . (3.5) 

Applying ( pJ.4D to ( ^.Ij ) we isolate the most general Gaussian part of the initial PI by 
performing an identical transformation as follows 

In particular, for the one-dimensional polaron we have 

aW[r] = aJdQoit, s, k) {e*^^^*'^) - e-^^'^^*"^)} - {T[r]) . (3.7) 

Since all Gaussian parts are concentrated, by definition, in the Gaussian measure, the new 



interaction ( |3.7| ) should not contain any quadratic path configurations. This requirement 



leads to the following equality 

= {T[r]) - ^JdQoit, s, k) e-^'^(*-) k' {[R{t, s)f) 

= i {{D-' - D,'),d) - aJdQoit, s, k) e-^^^(*-) k^ F{\t-s\). (3.8) 

Therefore, the adjustable function F{t) playing a key role in the new representation should 
be derived from the constraint equations: 



F{t)=- fdk [1 - cos(A;t)] D{k) , 
n J 



/ oo 

D{kU e + -^ fdt e-* [1 - cos(H)] F-^/\t) ] , (3.9) 

\ \/2ttJ I 







where D{k) is the Fourier transform of D{t). 

Exact analytic solutions to ( p.9|) are available in the weak- and strong-coupling limit: 



where 

A;2 7, (1 + ^)2 + \l-z\^ -2z^ 



1 2 °° 

7(A:) = 2v^ (Vi + v/TTP - v^) , xik) = ^jdz ■ 





Z3/2(1+;2)(P + (1 + ^)2) 



Note, function x(^) can be further simplified, but the result is a long expression consisting 
of rational and trigonometric functions. 

Substituting (|3^ ) into (|3.7| ) we rewrite the new interaction functional as follows 

W2[r] = JdQoit, S, k) e-'='^(*-^) hk[r{t)-r{s)]+k-F{t-s) _ ^ ^ :^[^(^) _ ^(^^pl (3^^^) 

Note again that W2 [t] does not contain any quadratic path configurations. Besides, the first 
cumulant is trivial (1^2^) = 0. 
Finally, we write 

^-mo^) = ( e«H^H ^^ ^ g-/3i^„{a) ^gaW.M^ ^ (3^^2) 

where -E'o(q;) is the general Gaussian contribution to the one- dimensional polaron GSE. 
Equations ( p.6D -( ^T2|) serve as the basis of the new representation which we call the 



generalized Gaussian approximation. The remaining non-Gaussian corrections should be 
evaluated by considering the following new PI 

J^(a) = (e"^2M) . (3.13) 



Note, by applying the Jensen- Peierls inequality to ( 3.13 ) one obtains 

J/3(a) >e"<'^^M> = l. (3.14) 

Therefore, Eo{a) represents a upper bound to the GSE of the one-dimensional polaron: 

Eo{a)>E{a). (3.15) 

ii. The described above scheme can be easily generalized to the case x 7^ that is 
necessary to evaluate the EM of the polaron. Omitting details of the extension we write the 
final result as follows 

6 



where the leading-order Gaussian approximation is: 



(3.16) 



g-/3£G(-,",/3X^xp J — [dk \n{k^D{k)) + - (d-' - D-\ D 
27iJ 2 ^ 



-(« dQ,o{t, s, k) e 



-k'^F(t-s) ikx{t-s)/l3 



(3.17) 



and the non-Gaussian correction is described by 



^aW2 lr,x] 



) = (exp < a dQ.o{t, s, k) 



-k'^F{t-s)+ikx{t-s)/l3 



^iklr{t)-r{s)]+k^F(t-s) _ ^ _^ — [r(t) - r(s) 



>)• 



iii. Now we consider the pure Gaussian approximation to Zf^{a,x) as follows 



Zp {a,x) = e 



-P£G{x,a,(3) 



(3.18) 



(3.19) 



Substituting (|3.17|) and (|3.19| ) into ( p.6|) we obtain the leading-order Gaussian approxima- 
tions to the GSE and EM as follows: 



Eo{a)=- — Idk 





In (^k'^D{k)) -k'^D{k) + l 

oo 

m:{a)=l + ^U fdtt'e-'F-'/\t) . 

2^/2^^ ^ ^ 



a 



'2t^ 



+ ^ dte-'F-'/\t) 



(3.20) 



Taking into account ( |3.10| ) we obtain the following analytic solutions 

. N _ / -a - (1/4 - 2/37r) a^ _ 0(^3) ^ a -^ , 

"^""^'l -ay 71-0(1), a^oo, 

f l + a/2 + {3/2-4:/7i)a'^ + 0{a^), a -* , 

1 (16/7r2)«4 + 0(a2), a^oo. 



mJa 



(3.21) 



iv. The Feynman estimate can be reproduced, if one builds a convex combination of the 
two known asymptotical solutions (|3.10| ) as follows 



Dp{k) = w/e + (1 - w)/{e + v^) 



(3.22) 



Substituting ( |3.22|) into (|3.2CI|) and optimizing the obtained energy with respect to pa- 
rameters {w,v} one is able to reproduce a Feynman-type upper bound Ep{a) for the 
one-dimensional polaron. It is inferior to the general Gaussian result for finite a, i.e. 
Epi^a) > Eo{a). Accordingly, m*p{a) deviates slightly from niQ^a). Note, however, that 
( p.22| ) is not the exact solution to 



IV. NON-GAUSSIAN CORRECTION 

The Gaussian leading-order terms Eo{a) and wi*(a) approximate well the exact GSE and 
EM of the polaron. To estimate the influence of the non-Gaussian correction we evaluate 
( p.l8| ) by using the following expansion: 

Jfs{a,x) = (e°^^['-'^l) = exp L {W2[r,x]) + y [{{W^ir^x])') - {W2[r,x]f] +..\ . (4.1) 

We stress that this is not a conventional perturbation series on the coupling constant a, 
each term of the exponent in the r.h.s. of ( [4.1| ) contains a in more complicated way and the 
result is a rapidly converging series even for large a. 

In doing so we restrict ourselves by estimating only up to the second cumulant in ( [4.1D . 
Appropriate analysis performed in the weak- and strong-coupling regimes for the third cu- 
mulant indicates that taking into account higher-order cumulants results in only tiny im- 
provement over the obtained estimate (see, e.g. ( |5.3|) and (|5.9| )). We suppose that this 
picture remains generally valid in the intermediate region of a. We obtain 

Jp{a, x) = exp I ~(3AE2{a) - ^AM2{a) - ... 1 , (4.2) 

where 

2/3 /3 oo 

AE2{a)=— lim — dtds dxdy dkdpp{t,s,x,y) , 

-oo 

AM2(a)=lim — dtds dxdy dkdp p{t, s,x,y) {k\t — s\ + p\x — yl)"^ . (4.3) 

-oo 

Above, we have introduced the following correlation functions: 

S(t, s, X, ^' F{t -x) + Fis - y) - Fis -x)- F{t - y) , 

p{t,s,x,i^ e-^'-'\-^''-y\ exp [-k'^F{t - s) - p'^F{x - y) + kpE{t, s,x,y)} . (4.4) 



The second-order non-Gaussian contributions in ( |4.3| ) may be derived analytically for the 
weak- and strong-coupling limit. For finite a we integrate out (|4.3| ) explicitly over variables 
k, p, X and y and the remaining double integrals are calculated numerically. 

Finally, taking into account both the leading-order Gaussian and the second-order non- 
Gaussian contribution we estimate the one- dimensional polaron GSE and EM as follows: 

E2ia) = Eo{a) + AE2{a) , 

m*2{a) = ml{a) + AM2{a) . (4.5) 



V. EXACT AND NUMERICAL RESULTS 

We have calculated analytically the GSE and EM of the one-dimensional polaron explic- 
itly with accuracy O(a^) and 0(1) in the weak- and strong-coupling limit, respectively 

Weak-coupling solutions 

In the weak-coupling limit, various perturbation methods give convergent series in powers 



of a. For example, analytic results up to the second order in a read jM 



Eps{a) = -a- (3/2v^ - l) a^ _ o{a^) , 

1 + a/2 + (5/8V2 - 1/A) a^ + 0{a^). (5.1) 



mps[a) 

A Lee-Low-Pines type method taking into account three-phonon correction to the Davydov 
phonon coherent state [O] results in the following numerical results 



4n 



Ecwwia) = -«-0.06066a2- 0.00844 a^ + 0(a 

rnhwwici) = 1 + a/2 + 0.19194 a^ - 0.06912 a^ + O(a^) . (5.2) 

By using the GGA method we calculate exactly the GSE and EM of the one-dimensional 
polaron up to the a^ terms. Knowing explicitly the weak-coupling solution ( |3.10| ) we calcu- 



late the second and third-order non-Gaussian corrections to the GSE and EM. Adding them 
to the leading-order Gaussian contributions we obtain 

Eo{a) + AE2{a) + AE^{a) 

= -a- (3/V8 - 1) a^ - (5 - 63^2/16 + 19\/3/16 - 29\/6/48) a^ - O(a^) , 
m*Q{a) + /S.m*2{a) + /S.m^{a.) 

= 1 + (1/2) a + (5/8V2 - 1/4) a^ + 0.0691096281 a^ + O(a^) . (5.3) 

Note that the remaining non-Gaussian corrections starting from the fourth order will modify 
only the neglected term O(a^). 

Strong- coupling solutions 

According to the Pekar "Produkt-Ansatz" ||^, the polaron ground-state wave function lip) 
in the strong-coupling regime {d = 3) is written as a direct product of the electron |\E') and 
field \(f) wave functions, besides, \(p) depends parametrically on |\E'). Further development 
of this method can be found, in particular, in [^,0 . A reliable numerical computation of 
the minimal value of the corresponding Hartree-type functional for d = 3 was performed 
in ||26|. The most rigorous investigations for the strong-coupling limit for d = 1 have been 



reported in [0,^,11]. Following the Pekar Ansatz the GSE of the one-dimensional polaron 
may be found by solving the following variational task 

F(a] T 

- lim ^.r^ = sup /c/x(^(x)A^(a;) + 2|^(x)rt . (5.4) 

^ ' ' —00 
The EM is obtained as follows: 



m*{a) 4 



oo 



a— >oo 



a"^ 71 



lim — ^ = - /rffcr |(^|e*''^|^)|^ (5.5) 



The corresponding optimized wave function obeying appropriate boundary conditions is the 
nontrivial solution of the following constrained differential equation: 

oo 

^"(a;)+4*3(a;)-^(a;) =0, /"rfx|^(a;)|^ = 1 . (5.6) 

— oo 

It admits an explicit analytic solution 

^(x) = V2/(e^ + e"^) = [72 cosh(a;)] "^ (5.7) 

that results in 

m*{a) = ^a^ + 0{l). (5.8) 

These solutions coincide (after the appropriate re-scaling) with the exact results reported 



first in 10 



Our leading-order Gaussian contributions in ( p. 21 ) differ from the adiabatic ones in (^]8| 



The reason is that as a increases the polaron self-interaction is less well approximated by a 
general Gaussian functional. Hence, non-Gaussian corrections are required to fill this gap. 
Taking into account higher-order corrections ( |4.1| ) we obtain a series converging rapidly to 
the exact value in (|5.8| ) as follows: 

E„{a) = -0.318310^2, 
Eoia) + AE2{a) = -0.327014 a^^ 
Eo{a) + AE2{a) + AE^ia) = -0.330205 a^ (5.9) 

Appropriate series for the EM reads 



<(«) = 


-- 1.621139 a\ 


m*(a) + AM2(a) = 


= 1.858065 a\ 


mlia) + AM^ia) + AMsia) = 


= 1.966430 a^ 



(5.10) 

Intermediate-coupling results 

For finite a we have solved equations (|3.9| ) numerically by means of an iteration accept- 
ing (|3.22D as the first approximation. After the sixth iteration step the numerical results 



do not change within the given accuracy. The obtained intermediate-coupling results for 
the Gaussian leading-order contributions and the second-order non-Gaussian corrections are 
presented in Tables 1, 2; and depicted in Figure 1 in comparison with several known data. 



10 



In conclusion, we have evaluated the ground-state energy and effective mass of the one- 
dimensional polaron by developing the generalized Gaussian approximation method. We 
have obtained exact analytic solutions to these quantities in the weak-coupling limit in- 
cluding the a^ term. For the strong coupling, we have derived the exact GSE and EM by 
adapting the Pekar adiabatic theory, then have estimated the same quantities within the 
GGA method by developing a systematic iterative scheme to approach rapidly the exact 
results. 

We have shown that for intermediate coupling the leading-order Gaussian contribution 
to the self-energy slightly improves the Feynman estimate and represents the lowest upper 
bound available. Considering the non-Gaussian corrections, we have calculated the second- 
order terms for the GSE and EM. Appropriate analysis performed for the weak- and strong- 
coupling regimes indicates that taking into account higher-order corrections results in tiny 
improvement over the obtained estimate. We deduct that this picture changes inconsiderably 
in the intermediate region of a. Our method works well in the entire range of a, does not 
require extensive numerical calculations and provides reliable results rather quickly. 

ACKNOWLEDGMENTS 

The author thanks G.V. Efimov for useful discussions. I am indebted to H. Leschke for 
reading the manuscript, helpful comments and for hospitality. 



11 



REFERENCES 

[1] L.D. Landau 1933 Phys. Z. Sowjetunion. 3 664 

[2] S.I. Pekar 1951 Research in Electron Theory of Crystals, (Moscow: Gostekhizdat); 1954 Un- 

tersuchungen ijher die Elektronentheorie der Kristalle (Berlin: Akademie-Verlag) 
[3] H. Frohlich 1954 Philos. Mag. Suppl. 3 325 
[4] R.P. Feynman 1955 Phys. Rev. 97 660 

[5] C.G. Kuper and G.D. Whitfield (eds.) 1963 Polarons and Excitons (London: Oliver and Boyd) 
[6] J.T. Devreese (ed.) 1972 Polarons in ionic crystals and polar semiconductors (Amsterdam: 

North-Holland) 
[7] Z.B. Su and L. Yu 1983 Phys. Rev. B27 5199 

[8] B.E. Sernehus 1987 Phys. Rev. B36 9509; 1988 Phys. Rev. B37 7909 
[9] A. Lorke, J. P. Kotthaus and K. Ploog 1990 Phys. Rev. Lett. 64 788 
[10] E.P. Gross 1976 Annals of Physics 99 1 
[11] H. Leschke and S. Wonneberger 1989 J. Phys.: Math. Gen. 22 L1009; Corrigendum 1990 J. 

Phys. A23 1425; 
[12] C. Quinghu, W. Kehn and W. Shaolong 1994 J. Phys.: Condens. Matter 6 6599 
[13] M.H. Degani and G.A. Farias 1990 Phys. Rev. B42 11950 
[14] F.M. Peeters and M.A. Smondyrev 1991 Phys. Rev. B43 4920 
[15] K. Okamoto and R. Abe 1975 J. Phys. Soc. Jpn. 31 1337 
[16] J. Poulter and V. Sa-yakanit 1993 in: Path Integration: Trieste 1991 eds. H.A. Cerdeira et al 

(Singapore: World Scientific) 506 
[17] J. Adamowski, B. Gerlach and H. Leschke 1980 in: Functional Integration, Theory and Ap- 
plication eds. J. -P. Antoine and E. Tirapegui (New York: Plenum) 291 
[18] M. Saitoh 1980 J. Phys. Soc. Jpn. 49 878 
[19] G.V. Efimov and G. Ganbold 1990 preprint DIAS-STP-90-22; 1991 NATO ASI series B255 

205 
[20] G.V. Efimov and G. Ganbold 1995 Phys. Part. Nucl 26 198 
[21] G. Ganbold and G.V. Efimov 1994 Phys. Rev. B50 3733 
[22] M. Dineykhan, G.V. Efimov, G. Ganbold and S.N. Nedelko 1995 Oscillator Representation in 

Quantum Physics (Lecture Notes in Physics, series Monographs m26 (Berlin: Springer) 1 
[23] G. Ganbold and G.V. Efimov 1998 J. Phys.: Condens. Matter 10 4845 
[24] N.N. Bogolyubov and S.V. Tyabhkov 1949 Zh. Eksp. Teor. Fiz. 19 256 
[25] S.V. Tyablikov 1951 Zh. Eksp. Teor. Fiz. 21 377 
[26] S. Miyake 1975 J. Phys. Soc. Jpn. 38 181; 1976 41 747 
[27] J. Adamowski, B. Gerlach and H. Leschke 1980 Phys. Lett. A79 249 



12 



FIGURES 

FIG. 1. The ground-state energy and effective mass of the one-dimensional polaron as func- 
tions of the electron-phonon coupling constant a. All curves are normalized to the corresponding 
Feynman estimates so that the horizontal dotted lines show the Feynman approximations. Dashed 
lines depict the leading-order Gaussian contributions and solid curves correspond to the corrected 
results taking into account the second-order non-Gaussian corrections. 
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TABLES 

TABLE L The ground-state energy of the one-dimensional polaron in the intermediate-couphng 
range. 



a 



0.5 
LO 
L5 
2.0 
2.5 
4.0 
6.0 
8.0 



En 



-0.5000 
-LOOOO 
-L5000 
-2.0086 
-2.7043 
-5.7934 
-12.155 
-21.067 



Ep 



-0.51006 
-1.04444 
-1.61314 
-2.23696 
-2.95968 
-6.04779 
-12.4074 
-21.3178 



Eo 



-0.51027 
-1.04532 
-1.61522 
-2.24044 
-2.96352 
-6.04964 
-12.4083 
-21.3183 



E2 



-0.51589 
-1.06710 
-1.66119 
-2.31150 
-3.05052 
-6.20592 
-12.7345 
-21.8871 



TABLE II. The effective mass of the one-dimensional polaron in the intermediate-coupling 
range. 



a 



0.5 
1.0 
1.5 
2.0 
2.5 
4.0 
6.0 
8.0 



m„ 



1.25000 
1.50000 
1.75000 
5.85933 
31.6793 
331.719 
1911.78 
6302.70 



m. 



1.32084 
1.88895 
3.11732 
6.83836 
21.4723 
281.622 
1784.90 
6068.77 



m„ 



1.32246 
1.89862 
3.15094 
6.90814 
21.4372 
281.474 
1784.82 
6068.73 



nio 



1.31091 
1.83258 
2.91361 
6.13458 
19.4943 
295.800 
1973.85 
6821.14 
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